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In this paper we investigate the relation between the scaling properties of the linear response func- 
tion R(t, s), of the thermoremanent magnetization (TRM) and of the zero field cooled magnetization 
. (ZFC) in the context of phase ordering kinetics. We explain why the retrival of the scaling properties 

' of R{t, s) from those of TRM and ZFC is not trivial. Preasymptotic contributions generate a long 

crossover in TRM, while ZFC is affected by a dangerous irrelevant variable. Lack of understanding 
of both these points has generated some confusion in the literature. The full picture relating the 
exponents of all the quantities involved is explicitely illustrated in the framework of the large A'^ 
model. Following this scheme, an assessment of the present status of numerical simulations for the 
Ising model can be made. We reach the conclusion that on the basis of the data available up to now, 
statements on the scaling properties of R{t, s) can be made from ZFC but not from TRM. From 
ZFC data for the Ising model with d = 2, 3, 4 we confirm the previously found linear dependence 
^ ' on dimensionality of the exponent a entering R{t,s) ~ s'^^^"^^ f{t/s). We also find evidence that a 

C/3 ^ recently derived form of the scaling function f{x), using local scale invariance arguments [M.Henkel, 

M.Pleimling, C.Godreche and J.M.Luck, Phys.Rev.Lett. 87, 265701 (2001)], does not hold for the 
I Ising model. 

S '. PACS: 05.70.Ln, 75.40.Gb, 05.40.-a 
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I. INTRODUCTION 



The behavior of systems out of equilibrium is a subject of wide current interest [1]. Most of the attention is focused 
on glassy or disordered systems. Nonetheless, many of the interesting features of slow relaxation, such as aging, can 
^ • be studied also in the simpler context of a phase ordering process. This is the dynamical process which takes place, 
' for instance, when a ferromagnet is suddenly cooled from above to below the critical point. Then, ordered regions 
\l grow by coarsening. The process is slow, i.e. the typical size of these regions grows with the power law L{t) t^/^, 
where z is the dynamic exponent. For dynamics with non conserved order parameter (NCOP), as it will be considered 
in this paper, z = 2 independent of dimensionality. In an infinite system equilibrium is never reached. Phase ordering 
• has been studied for a long time now [2]. However, despite its relative simplicity when compared to the complexity 
I of glassy behavior, still there remains lack of consensus and considerable confusion about the properties of the off 
"Y^ • equilibrium response function. This paper is devoted to clarify the issue. This is a problem not of minor importance, 

■ given that phase ordering is regarded as a paradigmatic example of out of equilibrium behavior. 
^ , For definiteness, let us think of an Ising ferromagnet with Hamiltonian 

^ : w = E (1) 

^ initially prepared at very high temperature and quenched at the time i = to a final temperature T < Tc- In a 
process of this type the initial magnetization is zero and remains zero at all times (ai{t)) = for < > 0. Quantities of 
interest are [3] the autocorrelation function 



C{t,s,to,tsc,teg) = {<J,{t)ads)) (2) 
where t > s > are two times after the quench and the linear (auto) response function 



R{t, S, to, tsc, te 



dhi{s) 



(3) 



where hi{s) is the external field conjugated to the order parameter. Traditionally, in phase ordering studies most 
of the attention has been devoted to the correlation function [2], while the response function has remained in the 
background. 
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In addition to the two observation times t and s, in Eq.s (1) and (2) we have exphcitly indicated also a dependence 
on the following characteristic times 



• to ^ A ^ . This is a microscopic time related, through the dynamic exponent z, to such a microscopic length as 
the lattice spacing or the inverse momentum cutoff A~^. 

• tsc- The process of phase ordering is characterized by dynamical scaling in the asymptotic time region (or late 

stage). The characteristic time t^c separates the preasymptotic from the asymptotic regime, i.e. it gives a 
measure of how much time is needed after the quench for scaling to set in. 

• tf^q. After the formation of domains of ordered regions, equilibrium is rapidly reached in the interior of domains. 
The characteristic time needed to establish this local equilibrium is the same as the equilibration time in the 
pure ordered phases. It is given by t^q ~ where ^ is the equilibrium correlation length in the pure phases at 
the final temperature T and z is the dynamic exponent introduced above. 

The correlation and response function can always be written as the s\im of two contributions [3] 

C{t,S,to,tsc,teq) = Cst{t — S,to,teq) + Cag{t, S,to,tsc) (4) 



R{t, S, to, tsc, teq) — Rstif ^ S, to, tpq) + Rag(t, S, to, t^c) (5) 

where the stationary contributions are what one has in equilibrium in the pure phases. Therefore, the usual fluctuation 
dissipation theorem is satisfied 

D /, J. J. \ ^ dCgtyt — S,to,teq) 

Rst{t- S,to,teq) = ^—g-^ -. (6) 

The rest, the aging contributions Cag{t, s, to, tsc) and Rag{t, s, to, tsc), are what is left over due to the existence of slow 
out of equilibrium degrees of freedom. The above split is useful for s suSiciently large, i.e. for 

s > teq (7) 

in order to have well separated time scales for equilibrium and non equilibrium behavior and for 

S > tsc (8) 

in order for Cag{t, s,to,tsc) and Rag(t, s,to.tsc) to exhibit scaling behavior. 
In connection with the aging contributions there are two basic questions 

i) how do Cag{t,s,tQ,tsc) and Rag{t,s,to,tsc) scale in the late stage 

ii) what is the relation between Cag{t, s, to, tsc) and Rag{t, s, to, tsc), if any. 

The second question belongs to the general area of the out of equilibrium generalization of the fluctuation dissipation 
theorem [4]. This is a problem not as trivial as it is believed to be for phase ordering systems [5,6], with interesting 
implications on the connection between statics and dynamics [7]. In this paper we concentrate on the first question 
which is preliminary to the second one. 

Assuming that s is large enough for (7) and (8) to be satisfied and dropping tgc, the scaling form of Cag{t,s,to) is 
given by 

Cag{t, s, to) ~ s-''g{t/s, to/ s). (9) 

It is well known [2] that 6 = 0. Furthermore, for s ^ to one can set y = in g{x, y) and it is also well known that for 
X 3> 1 one has g{x,0) = g{x) ~ x~^l^ , where A is the Fisher-Huse exponent. Information about Rag{t,s,to), instead, 
is scanty. Writing the scaling relation analogous to Eq. (9) in the form 

Rag{t, S, to) = to/s) (10) 
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both a and f{x,y) arc much less known than b and g{x,y). Despite considerable efforts, no consensus has been 
reached as of yet on the value of a. The situation for the scaling function /(x, y) is not much better. Recently, Henkel, 
Pleimling, Godreche and Luck (HPGL) [8], using local scale invariance [9], have derived an explicit form of the scaling 
function which is supposed to be of general validity. However, under close scrutiny this form appears neither to be 
obeyed in those cases where an exact solution is available, nor to fit numerical data for Ising systems, as it will be 
shown in section III. 

There is more than one reason for such an unsatisfactory state of affairs. The first one is due to a qualitative 
analysis [10] of the relation between the response function and the density of defects. A naive use of this argument 
leads to the conclusion that a is independent of dimensionality, e.g. for scalar systems a = 1/z. In this form, due 
to its simplicity, this argument has become deeply rooted in the literature [7,12,13], despite the accumulation of 
exact [14-16], approximated [5,13] and numerical results [5,17,18] incompatible with it. As we shall see, 5, to) is 

trivial in the sense that is proportional to the defect density only in the short time regime, but in no case this implies 
that a is independent of dimensionality. Another reason is that in simulations Rag{t, s,to) is too noisy to work with 
and, in order to deal with more manageable quantities, one must resort to the integrated response functions (IRF). 
The price for this is that reconstructing the scaling properties of Rag{t, s, to) from those of an IRF is not as simple as 
it might look at first sight [18]. This will be the main theme of the paper. 

We will show that, through the combined use of exact results and numerical simulations, definite conclusions can be 
reached for the exponent a by analysing in detail what actually goes on in the different methods employed to evaluate 
it. For what concerns the scaling function f{x,y), instead, our understanding of the problem remains incomplete. 

The paper is organized as follows. In section II we review existing information about Rag{t, s, to), we make general 
considerations on the scaling function and we comment on the HPGL theory. In section III we analyse the problem 
of retriving the properties of Rag{t, s, to) from those of an IRF concentrating on the zero field cooled magnetization. 
Section IV is devoted to the same problem from the side of the thcrmoremanent magnetization. In section V we use 
the solution of the large N model as an explicit illustration clarifying what goes on when different IRF are employed 
to obtain information on Rag{t, s, to)- Concluding remarks are made in section VI. 

II. WHAT IS KNOWN ABOUT Rag 

This paper is devoted to the study of the exponent a and the scahng fimction f{x,y) entering Eq. (10). We first 
summarise what is known from exact and approximate analytical results providing direct access to Rag{t,s,to)- We, 
then, make general considerations on f{x, y) and some remarks on the HPGL form for it. 

Ising model d=l 

In the exact analytical computation of the response function [14,15] in the d = 1 kinetic Ising model with Glauber 
dynamics, after taking s ^ tgc and neglecting to/s one finds 

Rag{t,s)^s'\t/S-1)-^''' (11) 

from which follows 

a = (12) 

and 

/(x,0)~(a;-l)-V2. (13) 
Furthermore, the correlation function is given by [19,20] 



C(M) = ^arcsiny'^ (14) 
which gives C{t, s) ~ {t/s)~^^^ for t/s^ 1. Hence, recalling z = 2, one has A = 1 and Eq. (13) can be rewritten as 
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It should be mentioned that a = has been found numerically also in the case of the kinetic Ising chain with Kawasaki 
dynamics [21]. 



Large N model 

Solving analytically the large N model we have found [16] (see also section V) Rag{t, s, to) of the form (10) with 

a={d-2)/2 (16) 

and 

„a+l-\/z _ 1 

where d is arbitrary and X = d/2. Notice that a = for d = 2. 
Gaussian auxiliary field ( GAP) approximation 

Berthier, Barrat and Kurchan [13] have calculated analytically an IRF using a GAF approximation based on the 
Ohta-Jasnow-Kawasaki method [2]. Prom their computation it is easy to extract Rag{t, s, to) which is in the form (10) 
with 

a=(d-l)/2 (18) 

and 

„a+l/2-A/z 

with A = d/2. Their calculation involves a diffusion constant of the form D = {d — l)/d which prevents letting 
rf — > 1, so they consider d > 2. Wo have worked out [5] an alternative GAF approximation, without restriction on 
dimensionality, which extends Eq.s (18) and (19) to d > 1. Then, we recover a = for d = 1 as in the Ising case. 

A. General form of f{x,y) and implications for Rag{t,s,to) 
All the above results for f{x,y) are of the form 

x->^ - e 

f{x,y)^- ; (20) 

where e = if the correlation length in the low temperature pure phase is finite, like in the d = 1 Ising model and in 
the GAF approximation, or e = 1 if the low temperature phase is critical [22] like in the large N model [16] (see also 
section V). 

We now make the phenomenological assumption that Eq. (20) is valid in general. Then the task becomes that of 
finding the exponents a, a and /3. For this it is useful to look at the short and long time behaviors. 

Short time behavior 

Let us rewrite Rag{t, s, to) introducing the time difference t = t — s in Eq. (20) 



Rag{t,S,to) = s"-^'+''^ 



{t/s + 1)-" 



(r + to)" 

Keeping r fixed and letting s to become large, to lowest order in t/s we find 



Rag{t,S,to) ~ S 



(r + to)" 



(21) 



(22) 
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with 

6={l + a)-{a-e) (23) 

and where e is the same as in Eq. (20). Therefore, from the short time behavior one can extract S. An important 
observation is that in the three exphcit cases considered above 6 coincides with the exponent entering in the time 
dependence of the density of defects. At the time s, this is given by 

L{s)-" ~ s""/^ (24) 

where L(s) and z are the domain size and the dynamic exponent introduced above, n = 1 for iV = 1, n = 2 for TV > 1 
and N is the number of components of the order parameter [2] . One can, then, immediately veriiy that 

S = n/z. (25) 

In the d = 1 Ising model and in the GAF approximation where e = and 

a = a + l/2 (26) 
from Eq. (23) we get 6 = 1/2, while in the large N model with e = 1 and 

a = a+l (27) 

we get 6 = 1. 

Long time behavior 

In the large time regime t/s ^ 1, from Eq.s (10) and (20) follows 

Rag{t,s,to)^s-^'+'^\t/s)-^'^/' (28) 

with 

Xr/z = a + l3. (29) 
Summarising, the exponents a, a and /3 can be obtained, in principle, by making three different measurements on 

Ragit, S, to): 

1. s dependence for fixed t/s gives a (from Eq. (10)) 

2. s dependence for fixed r gives S (from Eq. (22)) 

3. t dependence for fixed s gives Xr/z (from Eq. (28)). 

Before going into this, let us comment on the form of the scaling function derived by HPGL in Ref. [8]. 



B. Response function from local scale invariance 

Without making the separation (5) between stationary and aging components and neglecting the dependence on 
to, HPGL assume that the full response function R{t, s) obeys the scaling form 

i?(i,s)~s-(i+»)/HPGL(i/s)- (30) 
Then, using local scale invariance arguments they make the prediction that in general, for phase ordering, one has 

fHPGL{x) - _ ^^,+1 (31) 

where A is the Fisher-Husc exponent, provided there arc no long range correlations in the initial condition. In support 
of (31) they invoke the exact solution of the spherical model [23-25], which is equivalent to the large N model, and 
numerical simulations for the Ising model with d=2 and d=3. We make the following comments 
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• in the spherical or large N model Eq. (31) ideed reproduces the full response function. This coincides with 
the aging contribution (17) for x ^ I but not for a; ~ 1. This difference will turn out to be important (see 
section V). 

• Eq. (31) is contained in the general form (20) with e = 0, /3 = X/z — (a + 1) and a = a + 1. Inserting into 

Eq. (23) follows that in all cases one should have 5 = 0. Furthermore, one should also have a = a + 1 always, 
while from the explicit examples considered above this is true only in the large N case and not in the d = 1 
Ising model or in the GAF approximation. 

• HPGL theory is supposed to hold also for quenches at Tc- In that case the validity of Eq. (31) has been 
questioned in the framework of the field theoretic e-expansion for the response function [26] . 

• About the support to Eq. (31) from numerical simulations we will comment below. 

Now, in order to go beyond the explicitely solvable cases, the problem is to determine the exponents a, a and /3 in 
the Ising model with d > 1. As stated in the Introduction, measurement of Rag(t, s, to) is too noisy, so the program 
outlined above in subsection II A on the basis of Eq.s (22) and (28) requires an unrealistically long computing time. 
In the next section we discuss how to proceed with the help of IRF. 

III. ZERO FIELD COOLED MAGNETIZATION 

Indirect information on Rag{t,s,to) comes from numerical results on IRF. In general, an IRF is defined by 

t't2 

IJ,{t,t2,ti,to,tsc,teq) = / dsR{t, S,to,tsc,teq) (32) 
Jti 

with t > t2 > tl > and using (5) one has 

fl{t, t2,ti,tQ, tsc, teq) = fJ-stit — t2,t — ti,to, tgc, tgq) + IJ-agit, t2, tl,to, tgc)- (33) 

We will concentrate on the second contribution in the right hand side. The reason for introducing an IRF is that 

the integration over (ti,t2) lowers the noise. However, if one has to resort to an IRF, there is the related problem of 
retriving the properties of R{t, s) from it. This is not straightforward. IRF usually employed are 

1. the thermoremanent magnetization (TRM) 

p{t, tw,to, tsc, teq) = jJ-it, ^2 =twit\ = 0, to, tgcj teq) (34) 

obtained by looking at the response at the time t to an external field acting in the interval (0,f^) 

2. the zero field cooled magnetization (ZFC) 

X{t, tyj, to, tsc, teq) = IJ'it, t2 = t,ti = tw, to, tsc, teq) (35) 

obtained by looking at the response at the time t when the field acts in the interval {tyj, t). 

Both these quantities do have shortcomings. For TRM the problem is evident. The integration starts at ti = 0, so 
preasymptotic contributions are always included. The dependence on tsc cannot be neglected and this turns out to 
make it particularly hard to extract the asymptotic behavior in the cases of interest. 

With ZFC there is not such a problem. Taking t,^, ^ tgc, and neglecting tg,-, thereafter, one can be confident to be 
in the asymptotic region where scaling holds. So, using (10) with /(x, y) of the general form (20) and considering the 
case with e = 0, one has 

Xag{t, t^,to) = t-^Fit/t^M/t^) (36) 

with 

/•I ^l3+a-{a+l) 

F(x,y)=x-'' dz- -— . (37) 

^ ' ii/x (l-z + y/x)- ^ ' 



6 



The first observation is that if one seeks to determine a from Eq. (36) by looking at the behavior of Xag as is varied 
and X = t/tw is kept fixed, one must be aware that the dependence coming from to/t^ may play a role. In other 
words y = to /t^ may act as a dangerous irrelevant variable. Namely, defining the exponent a-)^ by 

Xag{t,t^) ^ t-'^-Xix) (38) 

there may be a difference between a in Eq. (36) and in Eq. (38). This depends on whether the integral in Eq. (37) 
diverges or not at the upper limit of integration as y — > 0. This, in turn, depends on the value of a. The second 
observation is that a can be extracted from the large x behavior of F{x, y), as we shall see in the following. Instead, 
the task of extracting f3 from Eq. (37) remains exceedingly complicated. 




A. The exponents a and 

The possibility that might not be identifiable with a, due to the presence of y = to/tyj, can be checked explicitely 
in the large N model [16] and in the GAF approximation [5,13], where ZFC can be calculated with arbitrary d. In 
both cases there is a value of the dimensionality such that y is dangerous irrelevant above d-^ . This implies that 
coincides with a for d < d-^ and is given by Eq.s (16) and (18). Instead, a-^ is different from a and is given by 

% = S (39) 

for d > d^ with S given by Eq. (25), which is independent of dimensionality. Logarithmic corrections appear Sit d = d-^, 
much in the same way as at the upper critical dimensionality in ordinary critical phenomena. The relation between 
a-)^ and a in these two models is given by 

for d < dy(^ 

with log corrections for d = d^ (40) 
for d> d-^ 

with 6 = 1 and d,^ = 4 in the large N model and with (5=1/2 and d,^ = 2 in the GAF approximation. We emphasize 
that in these two solvable cases, Eq. (39) holds only for d > d-^ where ^ a. 

Next, from extensive numerical simulations [17,5,18] of the Glauber-Ising model with d = 2, 3, 4 we have measured 
a-^ obtaining data which are fairly well consistent (Fig. 1) with the phenomenological formula 

(d-l)/4 for d<3 

1/2 with log corrections for d = 3 (41) 
1/2 for d>3. 

Since in the scalar case 6 = 1/2, it is evident that the pattern (40) is followed also in the Ising model with dy^ = 3. 

We may, then, conclude that in all cases: exact, approximate and numerical is given by Eq. (40) and that, 
therefore, the exponent a obeys the general formula 

where di is the dimensionality where a = 0. According to this picture, the distinction among the different systems 
comes through the values of (5, d^^, dL (see Table I). 





Ising 


GAF 


iV = oo 


s 


1/2 


1/2 


1 


dL 


1 


1 


2 


dx 


3 


2 


4 



TABLE I. 
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In this respect, notice that for N = 1 both from simulations and from GAF one has 6 = 1/2 and = 1, while 
there is a discrepancy between = 3 and = 2. However, this is not worrysome. As explained in Ref. [5], the 
dimensionality dependence of below takes place because d^ is the dimensionality below which minimization of 
magnetic energy competes effectively with minimization of surface tension in driving interface motion. Therefore, the 
balance of these two mechanisms is very sensitive to the treatment of surface tension and it should not come as a 
surprise that from an uncontrolled approximation, such as are those of the GAF type, a value of d^ which differs from 
the one observed in simulations is obtained. The shift from = 3 to d^ = 2 means that in the GAF approximation 
surface tension is overestimated with respect to simulations. 



0,75 



rt"^ 0,5- 



0,25 




FIG. 1. Exponent in the Ising model at various dimensionalities. The continous line represents Eq. (41), while the dots 
are the values from the exact solution of the model at d = 1 and from simulations with d = 2, 3,4. 



B. The scaling function x{^) and the exponent a 



Although the results described above [17,5,18] yield unequivocally d;^ = 3 for A'' = 1, in order to treat this point 
most carefully, we have investigated anew the behavior of Xagit^ tw) with very accurate simulations of the Ising model 
with NCOP, d = 2, 3, 4 and for different values of in order to get data also on the scaling function x(x), which has 
not been studied previously. 

First, let us illustrate the algorithm. There are several ways to isolate the aging contribution Xag{t,tw)- The most 
obvious is to compute the total xit^tw) by simulating a quenched system and then to subtract from it the stationary 
part Xst{t,ty^,) obtained by simulation of a system in equilibrium at the final temperature of the quench. A dilferent 
algorithm was introduced by Derrida [27] regarding the stationary contribution as due to thermal fluctuations inside 
the bulk of domains and the aging part as produced from the interfaces. The next step is to isolate the spins belonging 
to an interface. In order to do this a parallel simulation is performed of two systems with different initial conditions. 
The first is prepared in equilibrium at the initial temperature and then is quenched to the final temperature T, while 
the second is in equilibrium at the final temperature T from the beginning. These two systems evolve with the same 
thermal history at the temperature T. At each time step spins that are flipped in the first system but not in the 
second are considered as interfacial and their response is assigned to the aging part. 

These two methods arc equivalent, but also numerically very inefficient. Let us refer to them as global methods. 
The reason for inefficiency is that in order to extract the response produced by the spins on the interfaces one has to 
simulate the whole lattice. Since the interface density decreases as a huge amount of cpu time can be saved by 
an algorithm updating only the interfacial spins. We stress that a fast algorithm is crucial in order to have reliable 
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results in a numerically hard problem such as this. Therefore, we have adopted a no bulk flip algorithm, where a list 
of interfacial spins is updated at each move following the criterion that a spin belongs to an interface if at least one of 
the nearest neighbourghs is not aligned. Only moves of the interfacial spins are allowed. We then take the response 
of this system as Xag{t,tw)- 

In d = 1 it can be shown [14] that the no bulk flip algorithm corresponds to taking the limit of an infinite 
ferromagnetic coupling (J oo) in the Ising Hamiltonian and that this isolates exactly the aging part of the response 
function. With d> 1 the J ^ oo limit and the no bulk flip algorithm produce difl'erent dynamical evolutions and an 
argument analogous to the one in the d = 1 case cannot be made. What happens is that the limit J — > oo does not 
isolate Xag{titw) because, besides freezing spins in the bulk, it also freezes most of the interfacial spins. Notice that 
the no bulk flip dynamics does not obey detailed balance. This is simply due to the fact that bulk spins are frozen. 
However, this is not a serious problem since we already know that by restoring moves in the bulk detailed balance is 
recovered producing the stationary contribution in the response function, which we are not interested in. 

We have performed the simulations with the no bulk flip algorithm, after checking that the results are consistent 
within 5 per cent with those of the global algorithms. In practice, we measure the quantity 

1 

Xag {t,t^) = JjJ^ {(Ji ) hi (43) 

where hi is a quenched configuration of an uncorrelated random field, which takes the values ±/io with probability 
1/2. The angular brakets stand for the average over thermal histories, generated with the no bulk flip algorithm, and 
the overbar denotes the average over random field configurations. Simulations have been performed at T /Tc = 0.66 
for all values of d (for the lattice size and the number of realizations see Table II). Xag{t, tw) is measured in units 
and time in units of Monte Carlo steps. For each thermal history we have changed also the random field configuration. 
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1000 
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17000 










1750 


1024^ 


17000 










2500 


1024'' 


6000 











TABLE II. Lattice size Af and num- 
ber of realizations in the computation of 
Xag{t,tw) at different waiting times 



First, we have obtained by plotting Xag{t, tw) versus t^, for fixed values of x = t/tw In the range of explored 
there is excellent power law behavior. With a; = 7 we find (Fig. 2) a-^ = 0.28 for d = 2, = 0.47 for d = 3 and 
= 0.50 for d = 4. These numbers reproduce the results obtained previously [17,5,18] confirming that in the Ising 
model obeys closely Eq.s (40) and (42) with S = 1/2, d^ = 3 and di = 1 (see also Fig. 1). Furthermore, the observed 
behavior is with good accuracy independent of x, as it is shown in the inset of Fig. 2. The presence of a logarithmic 
correction at d = 3 is hard to establish from the data of Fig. 2 since we have only one decade in tyj. In Ref.s [17,5] 
where Xag{t,tw) was plotted against t for fixed tw over four decades, the logarithmic behavior is accessible. Also, it 
should be mentioned that Eq. (41) is a phenomenological formula, so it is hard to say whether the measured value 
= 0.47 for d = 3 is due to logarithmic corrections or to some other effect not captured by Eq. (41). In any case, 
the quality of the data for d = 2 allows to definitely rule out a-)^^ = 0.5, predicted by the qualitative argument referred 
to in the Introduction and to be discussed shortly. 
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10"^ 




B— H d=3 



10' 



10" 



10= 



10^ 



FIG. 2. Xasi'tj'tw) versus tw for fixed x — t/tw The slope yields a^. The plotted lines correspond to a; = 7. Values of for 
different values of x are depicted in the inset. 

Next, in order to investigate the scaling function x(a;) in Eq. (38), notice that from Eq.s (36) and (37) foUows the 
large x behavior 



X ° for a < 1 
x{x) ^ \ x^°'\ogx for a — \ 
x° 



(44) 



Using the values of a^^ from Fig. 2, we have plotted f^Xagit, tw) versus x = t/tu, for different values of (Fig.s 3,4,5). 
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FIG. 3. Scaling function £(a;) for the d = 2 Ising model with T/Tc = 0.66. 
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FIG. 5. Scaling function x(x) for the d = 4 Ising model with T/Tc = 0.66. 

Collapse of the data is obtained for x sufHciently large, where the scahng function decays with a power law and 
an exponent which coincides with a^. This is consistent with Eq. (44) only if a = a + 1/2 as in Eq. (26) and this 
rules out a = a + 1, which ought to apply according to the HPGL theory. Another way to see that a = a + 1 is 
untenable is that this would imply that x{x) goes to a constant for large x when a > 1. This, in turn, would lead to 
the unphisical conclusion that Xag{titw) does not decay to zero for large t and fixed t^] when d > d^. Therefore, we 
find that Eq. (26) holds for the Ising model not only for d = 1, but also at higher dimensionality. 

In conclusion, Eq.s (42) and (26) are our main results for a and a in the Ising model with d ranging from 1 to 4 
and with Eq. (40) explaining how a is related to a^. 
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C. Qualitative conjecture on 



We may now comment on the qualitative conjecture mentioned in the Introduction. Stating [10,7] that the aging 
contribution of ZFC ought to be proportional to the density of defects and assuming scaling, one finds 

Xag{t,tn,)^t-^X{tM (45) 

where S is given by Eq. (25). This requires ax = ^ ^ov all d contrary to the evidence presented above and summarised 
in Eq. (40), which restricts the validity of Eq. (39) to d > d^- This makes a big difference, for instance, in the d = 2 
Ising model where from Eq. (40) = 1/4, while from the above Eq.(45) follows = 1/2. In order to understand 
why Eq. (45) breaks down below d;^, let us go back to the behavior of Rag{t,s,to) in the short time regime. Prom 
Eq. (22) we may write 

Rag{t,S,to)^S-^h{T,to) (46) 

where /i(r, to) is some function of the time difference. The meaning of this is that the response, due to an impulsive 
perturbation at the time s, is proportional to the density of defects at that instant of time with a proportionality 
factor containing the retardation effect. This does not hold anymore in the long time regme r ^ s. When the time 
interval r is large with respect to s, multiple defect transits may have occurred through the observation site, spoiling 
the form (46). Sticking to the short time regime, i.e. taking t — tw <^tw and using Eq. (46), from the definition (35) 
follows 

Xag{t,ty,) t-\,{t - t^) (47) 

where Xs{t — tw) if^ a function of the time difference, which in Ref.s [5,17] we have identified with the ZFC associated 
to a single defect. Now, Eq.s (47) and (45) do require 

X{t/t^) ^Xs{t-t^) (48) 

which can hold only if both functions are constant. And this is precisely the point. As we have explained in Ref.s [5,17] 
Xsit~tw) contains the cumulative effect on a single defect of the perturbation acting all along the time interval (t^,, t). 
This saturates rapidly to a constant when the defect degrees of freedom act paramagnetically and the underlying defect 
motion is uncorrelated with the external field. However, at dimensionalities low enough to reduce surface tension below 
the threashold where the external field may take part in driving defect motion, Xs{t — i«,) acquires a non trivial time 
(iepenclence which renders S for d < d^- Finally, notice that in the framework of the qualitative conjecture with 

= S independent of dimensionality, there is no explanation for the exact d = 1 result = 0. Instead, according to 
Eq. (41) this exact result, far from being an anomaly, is embedded as a limiting behavior in the smooth dimensionality 
dependence for d < 3. 



IV. TRM 

Dealing with TRM, the separation (33) gives p{t, s) = Jq"" Rst{t — s) + Jg"' Rag{t, s). Contrary to what happens for 
ZFC, where Xst for long time saturates to a constant, here for the stationary contribution there are two possibilities: 
i) if Rst{t — s) decays exponentially Pst(t — tw) also decays exponentially or ii) if Rst(t — s) decays with a power law, 
like in the large N model, Pst{t — tyj) is subdominant with respect to pag{t,tyj). In both cases we can neglect pst and 
with it the distinction between p and Pag- 

As mentioned previously, TRM is affected by preasymptotic contributions which cannot be eliminated. This makes 
it quite difficult to establish if the asymptotic behavior has been reached in the simulations and ultimately to have a 
reliable estimate of a. In order to unravel what is the effect of the preasymptotic contributions on the scaling behavior 
of TRM, we have resorted as a guide to the solution of the large N model (section V). Here, we anticipate the results. 

Assuming tw > tsc, in the large N case there exists a dimensionality dp = 4 such that for d < dp TRM undergoes 
a crossover with a characteristic time t*, which may also be much larger than tsc- Introducing the effective exponent 



d\og p{t,tw,t* 



d log t 
one finds 



t/t ^^^^ 
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_ / for < t* , . 

""'-ff-Xa for (^0) 

For d = dp there is a crossover from a pure power law to a power law with logarithmic correction 



t^^'"E{t/U,) for <€. t* 

t~^'''\og{t^/tsc)E{t/t^) for t^»r. 



Finally, for d> dp one has the simple power law 

p{t,t^,t*)^t-^/'E{t/ty,) (52) 
and for all values of d, in the time regime considered the scaling function obeys 

E{x) ~ a;-^/^ (53) 

Taking this pattern as a guide (with dp, t* and exponents model dependent) let us now turn to simulations of the 
Ising model. Analysing data, the first thing to do is to check if a behavior of the type 

p{t,tj^t-''^p{t/tj (54) 

holds. If this is the case and if an exponent ap can be meaningfully extracted, the next problem is relating Up to a. 
According to the behavior found in the large N model, the identification Up = a can be made only if d < dp and 

t-w ^ t* . Numerical results for TRM in the Glauber-Ising model have been first obtained by HPGL [8]. Plotting 
p{t, tw) against x = t/tw for different tw in the range tyj G (25, 250) for d = 2 and tyj G (15, 100) for d = 3 they have 
obtained for Up a result of the form 

_ r 1/2 with log corrections for d = 2 , . 

'''' ~ 1 1/2 for d = 3 



and they have made the identification a = Up. 

The next round of simulations has been carried out by us [18] at the same temperatures and for the same system 
size as HPGL, but extending the range of tw up to 2500 for d = 2 and 250 for d = 3. Performing a different data 
analysis, i.e. plotting p{t,tw) versus for fixed x = t/tw, we have found good agreement between the slope of the 
curves in the large region, which in the log- log plot gives the cfFcctivc exponent (49), and the known values of A/z 
for the Ising model (A/z = 5/8 for d = 2 and \/z = 3/4 for d = 3). This is good evidence that in the scalar case TRM 
follows the crossover pattern obtained in the large N model when d < dp and with a crossover time t* larger than the 
we have reached in the simulations. Furthermore, on the basis of our data, we have estimated that 
the largest used by HPGL in Ref. [8] was not enough to enter the scaling regime (i.e. they had always tw < tsc) 
and therefore the values of ap they have obtained do not warrant any statement neither on the asymptotic value of ap 
nor on a. Our longer range of tw seems to be barely sufScicnt to enter the preasymptotic region where ap^eff = A/2;, 
suggesting that both d = 2 and d = 3 are smaller than dp, whose value in the Ising model, so far, we do not know. 
Hence, in order to observe the asymptotic exponent one shoud go to much longer waiting times tw- 

Henkel and Pleimling [28] have produced new simulations for d = 2 extending the range of tw up to 5000. Plotting 
p{t,tw) versus tw for fixed x and adhering to the point of view that the TRM data are aff'ected by a long crossover, 
they claim i) to have succedcd in going past the crossover time reaching the asymptotic region and ii) to have found 
that Eq. (55) is verified. The objection to this claim is that in d = 2 one has \/z = 5/8 > 1/2 > a^ = 1/4. Therefore, 
even if a decrease of the slope from a number close to X/z = 0.625 toward 0.5 is observed over a narrow time window, 
there is no way to decide whether the true asymptotic regime has been reached or the slope might still keep on 
decreasing, by going further with tw, untill reaching asymptotics at 0.25. 

In other words, the new simulations in Ref. [28] leave the issue undecided and yet longer simulations are needed. 
Despite, by now, there is sufiicient evidence that TRM is not the most efiicient and reliable way to get to the exponent 
a, we have undertaken new simulations with tw up to 5000 for d = 2 and 500 for d = 3 (for the lattice size and number 
of realizations see Table III). 
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realiz 
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25 


1024^ 


2000 


100^ 


1500 


50 


1024^ 


2000 


150^ 


2500 


100 


1024^ 


2000 


150^ 


2500 


ZOU 


1024^ 


13000 


150^ 


2500 


500 


1024^ 


16000 


160^ 


2500 


1000 


1024^ 


18000 






1750 


1024^ 


23000 






2500 


1024^ 


13000 






5000 


2048^ 


7000 







TABLE III. Lattice size JV and num- 



ber of realizations in the computation of 
pag{t,tm) at different waiting times 



The double logarithmic plot of p{t,tw) versus for fixed x shows (Fig. 6) that a power law behavior, possibly, 
sets in only in the region of the largest tw reached. Taking the slope in this region as a measure of a^ff^p we find 
values (inset of Fig. 6) which lie above 0.5 for all x and that are just below A/z = 0.625 for d = 2 and A/z = 0.75 for 
d = 'i. Hence, although we have reached the same maximum value of as in Ref. [28] for d = 2 and we have gone 
much farther for d = 3, we may state that no evidence of asymptotic behavior with Op = 1/2 is found. Rather, the 
combination of the d — 2 and d — 3 data in Fig. 6 shows unequivocally that, at best, only the onset of the scaling 
region is entered where ap,eff is about to take the preasymptotic value A/z, confirming the picture obtained in our 
previous work [18]. 




d=3 



■^3 4 

10" 10 10 

t 

w 

7. The slopes in the large region yielding ap^^ff for different values of x are 

In summary, we have accumulated sufficient numerical evidence to establish that TRM data fit in the general 
pattern of behavior abtained from the solution of the large N model, with dp > 3 and a value of t* which is greater 
than the largest t^, reached so far. Therefore, since asymptotics has not been reached, no statement on a can be made 
from the present knowledge of TRM. 
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FIG. 6. pacj{t,tw) versus for fixed x 
depicted in the inset. 
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Finally, let us make a comment on the quotation in Ref.s [8,29,30] of the analytical solution of the GAF approx- 
imation by Berthier et al. [13] as a support to the claim that a is given by Eq. (55). In fact, here is where is most 
evident the type of confusion that can be made by not being careful about which exponent one is talking about. In 
their computation Berthier et al. find = 1/2 for d > 2 with logarithmic correction at d = 2, as in Eq. (55)which, 
however, is meant for a. What one should have clear in mind is that they compute an for d> d^, i.e. right where 
7^ a. This can be checked recalling that in the GAF approximation a is given by Eq. (18) and that d^ = 2. Hence, 
for d = 2 the logarithmic correction belongs to and not to a. For d = 3 it is that takes the value 1/2, while 
from Eq. (18) follows a = 1. So, the results of Berthier et al. certainly cannot be quoted if one wants to identify with 
a an exponent obeying Eq. (55). 



V. TRM AND ZFC IN THE LARGE N MODEL 

In this section we study in detail the large A'' model [16,31] as a useful example which gives the complete picture of 

what happens when looking at the different response functions introduced above. 

Consider a system with vector order parameter (j){x) = {4>i{x), (j)N{x)) and Hamiltonian of the Ginzburg-Landau 
form 



H[0] = j d'^2 



(56) 



where r < 0, g > 0. In the large N limit the equation of motion for the generic component of the order parameter in 
Fourier space is given by [16] 

^^^M = -[k^ + I^tmk^t)+v{k,t) (57) 
where ri{k, t) is a gaussian white noise with expectations 

{■n{k,t)r]{k',t')) =2T{2T^Y6{k + k')5{t-t') ^ ' 

T is the temperature of the quench and the function of time 

I{t)=r + g{^^ix,t)) (59) 

must be determined self-consistently, with the average on the right hand side taken both over thermal noise and initial 
conditions. The formal solution of (57) is given by 

(pik, t) = R{k, t, O)0o(fc) + / dsR{k, t, s)r?(fc, s) (60) 

Jo 



where R{k, t, s) is the response function 



R{k,t,s) = ^e-'''(*-^^ (61) 



with Y{t) = exp[Q{t)] and Q{t) — dsl{s). With an uncorrelated initial state at very high temperature the initial 
condition 4>o{k) = 4>{k,t = 0) can be taken to be gaussianly distributed with expectations 



km = 

h{k)Mk')) = M^^)''S{k + k') 



(62) 



The actual solution is obtained once the function Y{t) is determined. In order to do this notice that from the definition 
of Y{t) follows 



dV^t) 
dt 



= 2[r + g{<i>\x,t)]Y\t). (63) 
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Writing {(jr{x,t)) in terms of the structure factor 



(0^(^,i))= I ^,C{k,t)e-''/^' (64) 



where A is the momentum cutoff and using Eq.(60) to evaluate {(j){k,t)(j){k' ,t)) = C{k,t){2iT)'''6{k + k') we obtain 



C{k,t) = R'^{k,t,0)A + 2T [ dsE?{k,t,s). 

Jo 

Then, inserting Eq. (64) into Eq. (63) we obtain the integro-differential equation 

^ = + 2,A/ (t + ^) + 4,r| <<./('- + 5^.) y'is) 

where f{x) = J ^^J^d ^~'^^^^ = (Sttx)"^. After solving this equation, the response function is given by 



(65) 



(66) 



R{t,sM) = J ^i?(fc,i,s)e-'=VA= = (4;r)-'^/2^(t-s + to)-'^/' (67) 

where to = 1/(2A'^). 

Let us now come to the identification of the general structure of Eq. (5). Since in the stationary regime Y{t) is 
time independent, we immediately obtain 

Rst{t - s,to) = (4^)-^'/2(t - s + to)-''/' (68) 

and 

' (i-s + io)-''/'. (69) 



Rag{t,S,to) = {4Tr)-''/^ 



Y{s) _ ^ 



Yit) 



Notice that Rsi (t — s, to) is temperature independent, implying that there is a stationary response also at T = 0. This 
holds for soft spins, while for Ising spins there is no stationary response at T = 0. 

Next, in order to investigate the scaling properties we must first learn about the time dependence of Y{t). We do 
this in the T = case, since quenches below the critical point are controlled by the T = fixed point [2] . Making the 
ansatz Y{t) = At~'^ from (66) one gets 

Au;t-(2-+i) = rAt-'^ + j^^{t + to)-"/^ (70) 

and assuming that the left hand side is negligible one finds w = ci/4 with A = {8it)~'^^'^A/Mq, where Mq = \J —rjg 
is the zero temperature magnetization. This is consistent if, in addition to t ;» to, one has also 

t > t,e = -rf/(4r) (71) 

where the characteristic time tgc sets the time scale over which the three terms in Eq.(70) are all of the same order 
of magnitude. Therefore, t^c is the characteristic time separating the early from the late stage. 

The above described behavior of y(t) is illustrated in Fig. 7 displaying the numerical solution of Eq. (66) for 
different values of r. In all numerical computations we will take A = 1, T = and time is measured in units to. The 
onset of the scaling behavior is sharp and we have identified tgc with the time where the power law begins (inset of 
Fig. 7). Then, for s,t > tgc from Eq.s (67) and (69) we have 

i?(t,s,to) = s-(i+'^)7(t/s,to/s) (72) 

with 

]{x, y) = {Awy/^x^ix - 1 + y)-(i+«) (73) 

where 
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a=(d-2)/2. (74) 

as in Eq. (16). 

The connection between lo and A/z can be established from the autocorrelation function. Keeping on considering 
T = 0, from C(fc, t, s) = R(k, t, 0)R{k, s, 0) A follows 

C{t, to) - / ^Cik, t, s)e~^''^' (75) 

= (47r)-''/2^s2"-(i+")(t/s)'^[t/s + 1 + to/s]"^^+°^. (76) 
The requirement limt^oo G{t, t) = Mq implies 

2w = 1 + a (77) 

and comparing Eq. (76) with Eq. (5) we find 

uj = X/z. (78) 
Hence, in the large TV model, A and a are not independent exponents, since from Eq.s (77) and (78) follows 

A=l + a. (79) 
Nonetheless, for generality we shall keep on using the notation with two different exponents A and a. 




FIG. 7. Y{t) for difTcrent values of r and T — 0. tsc is estimated at the onset of the power law behavior and plotted against 
r in the inset. 

Finally, for the aging contribution (69) we may write 

Ragit, S, to) = ^o/s) (80) 

with 

fix, y) - (47r)-'^/2(x" - l)(x - 1 + j/)-(i+'^' (81) 

and writing lu = 1 + a — X/ z Eq. (17) is recovered. 

The above result shows that in the large N model it is not only Rag{t, s) to scale, but also the full autoresponse 
function R{t, s). This, obviously, means that i?st(t — s) obeys scaling, as it can be checked immediately from Eq. (68) 
and this is a consequence of the fact that the whole low temperature phase is critical. 
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A. TRM 



We now explore the properties of the IRF in the large N model. Let us begin from TRM. Since the explicit 

forms (68) and (80) with (81) show that Rst{t — s, to) decays faster than Rag{t, s, to) with the time separation t — ,s, 
taking t ^ tyj and using the definitions (33) and (34) the stationary contribution to TRM can be neglected. Hence, 
in the following we will ignore the distinction between p and pag- Furthermore, taking tyj > tgc and dropping the 
dependence on to we can write 



P(t, tyj^ tgc) 



dsY{s) 



dsY{s) 



(82) 



where wc have separated the preasymptotic from the asymptotic contribution in the integral. Wc shall sec shortly 
that the first one plays a crucial role. Introducing the notation B{tsc) = Jq"" dsY{s) and using Y{s) = ^g-^/^^-Ca+i) 
in the second integral, we find 



P(t^ tyj^ tgc) ^1 



-A/z 



Ko + K^tl'^-'' 



where 



(47r) 



-d/2 



B{tso) 



t 



X/z—a 



{\/z - a) 



and 



{A-nY/-^{X/ z - a)' 



(83) 



(84) 



(85) 



Eq. (83) is the main result from which follows the non trivial dependence of Op on dimensionality. Notice, that all the 
dependence on the preasymptotic behavior is collected in Kq and the very presence of this non vanishing term entails 
that the asymptotic power governing TRM is either X/ z or a according to the sign of {\/z — a). Therefore, writing 
\/ z — a = {dp — d) / dp with dp = 4 we have a crossover for d < dp, logarithmic corrections for d = dp and a correction 
to scaling for d> dp. 

Introducing the characteristic time 



\Ki\ 



Eq. (83) for d ^ dp can be rewritten as 



p{t,t^,t*) = t-''''E{t/t^,t* /t^) 



with 



E{x,y)=Ko [l±y^/^-» 



-A/2 



(86) 



(87) 



(88) 



where the + and — signs apply to d < dp and d> dp, respectively. In the first case the crossover time t* is given by 
(Appendix) 



(89) 



showing that t* is a new time scale which can become much larger than tgc- Instead, in the second case from Eq. (86) 
follows (Appendix) 



t*/tsc < 1 

implying tyj/t* > 1 for any tyj > tgc- Finally, for d = dp from Eq. (83) we have 



p{t,t^,t\tsc) = [A-Ky/'t. 



d/2,-\/z 



1 + 



\0g{t*/tsc) 



l0g{t*/tsc){t/t 

w J 



-X/z 



(90) 



(91) 



18 



where t* is given by 

t*/t,c = e^*- (92) 

and C is a constant (Appendix). 

Therefore, as anticipated in section IV, the scaUng properties of TRM exhibit the following dimensionality depen- 
dence 



d < d„ 



There is a crossover with the effective exponent 

dlogp{t,t,j„,t*) 



= X/z 



91ogt„ 



1 + (tw/t*)^/'-" 



(X/z-a) 



yielding 



_ X/z for tyj < t* 



(93) 



(94) 



The crossover involves a logarithmic correction 



t^u^^^E{t/t^) for < t* 

t~^''\og{t^/tsc)E{t/t^) for t^»r. 



d > dp 

There is a pure power law for all > isc 

p{t,t^,t*)=t-^''E{t/U) (96) 

with 

E{x)r^x-^/\ (97) 
In the end, in the large A'' model the relation between a and the exponent ap appearing in Eq. (54) is given by 

{a for d < dp 
X/z with log corrections for d = dp (98) 
X/z for d> dp 

where dp = 4. 

In order to illustrate the behavior of TRM we have solved numerically for p{t,tyj). In Fig. 8 we have plotted the 

effective exponent (93) versus for different values of r (giving rise to different values of tsc), with fixed x = t/tw = 20 
and for d = 2.1 < dp. The curves show quite clearly three different regimes: the early regime to the left of the peak 

— X/z 

followed by the intermediate regime going like tw , whose size depends on tgc, and eventually by the late stage 
regime going like t~°' ■ The value t^ax of at the peak can be identified with tsc since it depends on r according to 
Eq. (71) (see inset of Fig. 8). For completeness we have plotted the same figure for d = 5 > dp (Fig. 9) which shows 
the existence only of the early regime followed immediately by the asymptotic regime with the exponent X/z (without 
any crossover or intermediate scaling regime) according to Eq. (96). 
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FIG. 9. Effective exponent a^.e// in the large A'^ model versus for different values of r with x = 20, d = 5 and T = 0. 



B. ZFC 

Taking > tsc and using the definitions (33), (68), (69) we have 
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2^1-<i/2 

Xstit - t^,to) = (4^)d/2(^_2) {1 - - + 1]""} (99) 

and 

Xag{t, U,to) = t'Tit/t^, to/t^) (100) 

with 



F{x,y) = {Awy/^x-" f (iw(tt^/"-(i+") - 1)(1 

Jllx 



y/x)-^'+''\ (101) 



Therefore, in order to estabUsh how Xag scales with f^, it is necessary to know how the scahng function F{x, y) behaves 
for small y. As already pointed out, this depends on the behavior of the integral at the upper limit of integration, 
which is convergent (divergent) for a < 1 (a > 1). Hence, from 1 — a = {d^ — d)/2 with d-^ = A follows 

.T"" for d < 

F{x,y) ^ { x-^logix/y) for d ^ d^ (102) 
y^~"/x for d> d^. 

Inserting into Eq. (100) and comparing with Eq. (38) we recover Eq.s (40) and (42). Finally, for large x we obtain 
the analogous of Eq. (44) 

x~°' for d < d^ 
X{x) ^ { x-^logx for d = dy^ (103) 
for d > d^. 

Notice that the separation of the stationary from the aging response function has played a crucial role. Had we used 
the form (31) of HPGL in Eq. (101) we would have obtained a completely different behavior, with d^ = 2 and in place 
of Eq. (40) 

a = rf/2 - 1 for d<2 
a^ = ^ with log corrections for d = 2 (104) 

for d> 2. 

In order to illustrate the difference in the behaviors of TRM and ZFC we have solved numerically for Xag{t, tw) and 
for the corresponding effective exponent a^^ef fitw, (Fig-s 10,11) with the same values of d and r used for TRM. 
These Figures show that both above and below d^ there is no crossover, but there is only the early regime followed 
abruptly by the asymptotic power law behavior, as for TRM above dp (Fig. 9). Furthermore, we have depicted in 
Fig.s 12,13 the scaling function x(x), obtained by plotting t^Xagitw,^:) versus x for different t^, which obeys the 
power laws (103) for large x. These arc the analogous of Fig.s 3,5. 

We can now summarise what we have learned from the large N model about the connection between a^, a^^ and a. 
In this case the explicit solution (74) is available and a is a linearly increasing function of dimensionality vanishing 
at dL = 2. The question is how much of this could have been inferred relying only on the information from TRM or 
ZFC. The answer is that both Op and coincide with a below certain dimensionalities dp and d^^. At d — dp and 
d = d^ there are logarithmic corrections, while above these dimensionalities ap and are different from a and differ 
one from the other (Fig. 14). Although in the large N model dp = d^ = 4, we have kept distinct notations because 
dp, which is the dimensionality where X/z = a, and d^ the dimensionality where a — 1 = 0, need not to coincide 
in general. In the large N model they do coincide because of Eq. (79). Furthermore, even below dp and d-^, where 
Up = a^ = a, there remains a considerable difference between TRM and ZFC in relation to the time scales {t* and tsc) 
over which these exponents are observable. Comparing Fig.s 8,10 one can see at glance that the difference between 
these time scales in certain conditions, here set by the value of r, can become huge and if working with TRM it may 
require an enormous tu, before reaching the asymptotic regime where Up and a can be identified. 
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FIG. 10. Effective exponent a^.e/f in the large A'' model versus tm for different values of r with x = 20, d = 2.1 and T = 0. 
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FIG. 11. Effective exponent a^.eff in the large A'^ model versus for different values of r with x = 20, d = 5 and T — 0. 



22 




FIG. 12. Scaling function x{^) in the large A'' model with d = 2.1 and T — 0. In this case ax = a = 0.05. 
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VI. CONCLUDING REMARKS AND OPEN PROBLEMS 



In conclusion, we have shown that all existing analytical results and the numerical evidence coming from ZFC in 

the Ising model are consistent with an exponent a of the form (42). The dimensionahty independent behavior (39) 
predicted by the qualitative argument for holds only for d > where ^ a. This is due to the presence of 
a dangerous irrelevant variable. Once this is taken into account, analytical and numerical results form a coherent 
picture and the issue can be considered as settled. 

For what concerns Eq. (55), regarded in Ref.s [8,28-30] as the exponent a in the Ising model, we have shown that it 
does not have any analytical foundation, because Ref. [13] contains a computation of a^. Furthermore, the numerical 
evidence, being based on TRM data, is inconclusive since the largest tw reached so far are below the crossover time 
t* . Therefore, t„ is still far from being well inside the asymptotic region as required for the TRM data to qualify as a 
challange to those obtained from ZFC. There is no doubt that among all possible IRF that one can employ to study 
the exponent a, TRM is the most unfavourable and the less reliable one, as abundantly explained in the previous 
sections. 

For what concerns the scaling function f{x, y), our ZFC data are consistent with an f{x, y) in the Ising model of the 
form (20) with the exponent a = a + 1/2 in place of a = a + 1, appearing in the HPGL theory. We have also shown 
that with the HPGL theory it is not possible to reproduce the short time behavior of Rag{t,s,to). Nonetheless, our 
knowledge of the scaling function f{x, y) is still incomplete, since from ZFC data we cannot determine the exponent 



After this survey of what can and what cannot be done with ZFC and TRM, it seems clear that in order to study 
Ragit, s, to) the right thing to do would be to use neither of them. Rather, one should use an IRF of the general form 
(32) with ti >• tsc to eliminate the crossover affecting TRM and with t2 <t'm order to avoid the dangerous irrelevant 
variable in ZFC. Namely, assuming the form (20) of f{x, y) and using Eq. (32) one should consider 



If t2 <t and to, the dependence on to can be neglected and the above equation can be used in two ways. Rewriting 



and keeping xi = ti/t and X2 = t2/t fixed, the exponent a can be measured. Next, for t^t2 from Eq. (105) follows 




(105) 




(106) 
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Jti/t 



(107) 



from which /3 + a can be measured, while a as we have seen can be extracted from x(a;). We plan to pursue the 
investigation of this IRF in future work. 

Finally, the results obtained in this paper open a number of interesting problems in the general theory of phase 
ordering. We stress that our results are phenomenological. In particular, we do not know why dL and take the 
values they take, rf/, seems to coincide with the ordinary lower critical dimensionality, but we do not know whether 
this is really so, or it is just a coincidence. Even less we can tell about the values taken by the upper dimensionality 
d^. It should be noted the failure of the GAF approximation to reproduce the correct dependence of a on in the 
scalar case. In short, we have no theory for the observed behavior of the response function in phase ordering kinetics. 
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VII. APPENDIX 



From the definitions (84), (85) and (86) we have 

^ \B{tsc)\X/z-a\ 



l/{X/z-a) 



(108) 



with the +(— ) sign if X/z — a < (>)0, i.e if d > (<)4. In order to estimate B{tsc) we use the linear approximation 

ftsc 



B{tsc)= [ 

Jo 



[e^^- - l)/r. 



Then, using A = {Sn)-'^''^ A/M^ and t^c = -d/Ar we find 



with 



and inserting into Eq. (108) we get 



B{tsc)/A = Ctsc 
C = 4(87r)'^/2(l - e-'^/^)M^/Ad 



t*lts 



(4-d) 



ctit" ± 1 



4/(4-d) 



(109) 



(110) 



(111) 



(112) 



For d < 4 the above equation must be taken with the minus sign. This requires tg^ > [C(4 — (i)/4]~^/'^ or 4|r|/d < 
[C(4 — d)/4]^/'*. To lift this restriction on the value of r one must do better than the linear approximation in the 
estimate of B{tsc)- Taking tgc large enough Eq. (89) is obtained. 

If d > 4, instead, from Eq. (112) follows t* < tgc justifying Eq. (96). 

Finally, for d = 4 from Eq. (83) we get 



pit-, tyj^ t , tsc) 



-XI z 



(47r)-''/2(t/t^)-V^ 



and defining t* by 

\og{t*/U,) = B{Uc)/A 
Eq.s (91) and (92) are recovered after using Eq. (110). 



(113) 



(114) 
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